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Abstract 

The problem of the p meson resonance solution in the extended NJL model is 
studied. It has been shown previously that the solution for the p meson changes 
smoothly from a bound state to a "virtual state" when the vector channel attraction 
becomes too weak to support bound state solutions. We have found that together 
with these solutions a separate resonance branch appears which is not connected 
smoothly with the bound state solutions. At small values of constituent quark mass, 
m ~ 200 MeV, it corresponds to the physical vector meson resonance solution in 
the quark-continuum region. The a\ and a resonances are also discussed. 
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1. We have become accustomed to think of resonance states in the extended Nambu- 
Jona-Lasinio (NJL) model [IJ as being solutions of the corresponding Bethe-Salpeter (BS) 
equations, analytically continued into the complex p 2 -plane. This follows from our quan- 
tum mechanical experience in describing a system which can disintegrate. Such a system 
can have a very small disintegration probability and is then known as a quasi-stationary 
or a resonance state. The propagator of a quasi-stationary state has always poles in the 
complex po-plane which are situated close to the real axis. This method can also be ap- 
plied to systems with a relatively large disintegration probability. Conceivably, we deal 
with such a kind of states in the case of the strongly interacting quark- ant iquark systems 
in the extended NJL model, if they are located in the continuum region ||. 

The question of the location of a complex pole which is close to the real axis can 
actually be resolved without going in the complex p 2 plane. In this letter we show how 
to do that. This simple procedure is quite general and can be used whenever one is 
interested in narrow resonance solutions of a known propagator or self-energy function of 
the corresponding state. The idea to find out the resonance parameters without going 
into the complex p 2 plane has been used in the spectral-function (SF) method It 
has been assumed that a maximum in the imaginary part of a meson propagator can be 
associated with a resonance mass. This approach has been criticized in || . It was argued 
there that although the SF method and the analytical continuation of the BS equations 
into the complex p 2 plane yield similar results for the a\ meson, the same methods lead 
to qualitatively different descriptions in the case of the p meson. Using our approach it is 
not difficult to understand why the SF method under certain circumstance may lead to 
wrong conclusions. 

There are several reasons why the resonance solutions in a quark continuum region are 
important. First, the a\ meson always lies above the qq threshold. Second, the a meson is 
displaced in that region when the current quark mass is different from zero. Furthermore, 
the statement that there is no resonance solution for the p meson [fj, |6|] requires further 
investigation. In fact, it can not be excluded that for solutions of the relativistic BS 
equation there exists a similar difference in the behaviour of s- and higher partial waves 
as it is the case in nonrelativistic potential theory. However, it is conceivable that the 
relativistic BS equation exhibits some new features in comparison with a nonrelativistic 
theory. In this letter we will show that besides the bound and "virtual state" || solutions, 
which are smoothly connected at some value of the vector interaction strength Gy, there 
exists for all values of Gy a resonance branch for the p-meson. We show that an overall 
fit to empirical data can be achieved either with large constituent quark masses (~ 400 
MeV) or small ones (~ 200 MeV). In the former case the physical p meson appears as 
a bound qq state. The resonance solution corresponds to the physical p meson for small 
constituent quark masses. This finding lends credit to the approaches which use small 
values of non-strange quark masses m, i.e. when m p > 2m, as it is frequently used in the 
literature ||. It also shows that even so in the NJL model confinement is absent, one can 
still investigate some properties of states which a priori seem to lie outside the range of 
applicability of the model. 
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2. Let the propagator of some state be defined by 



A( P 2 ) = rV) • (i) 

It frequently happens that the function f(p 2 ) becomes complex-valued for certain values 
of external momenta, and it is thus conceivable that the pole will be displaced into the 
complex p -pla,ne. One assumes usually that near the pole the propagator behaves as 

where 7 is positive, so that for infinitesimal small 7 we make contact with the standard 
ie prescription for the propagator poles. In this case the complex pole coordinates p 2 = 
M 2 — 27 are interpreted as a mass, M, and a width, T = j/u (p), of a corresponding 
resonance state. Here u(p) = \/{p 2 + M 2 ) is an energy. This interpretation follows from 
the fact that the leading contribution of the Fourier-transformed two-point function at 
large time originates from the momenta associated with the propagator pole in the lower 
complex po-plane. 

If 7 is small (7 <C M 2 ) the equation f(M 2 — ij) = can be solved in the small width 
approximation, that is by expanding it around the point M 2 — iO in a Taylor series and 
restricting oneself to the terms linear in 7. This is possible if the function / is analytic 
everywhere inside a circle Co with center at M 2 — iO and M 2 — 27 G C . We assume that 
the branch cut of / does not pass the circle Co- In this case one obtains 

7=4", R f R' f + I f I' f = 0, (3) 



R' f 



where 



R f = Ref(p 2 ), I f = lmf(p 2 ), 
R , = dRefjp 2 ) ^ r dlmfjp 2 ) ^ 

taken at p 2 = M 2 . Solving the second equation in (§) one can obtain the mass and the 
first one gives the width. It is, however, mandatory to check that the self-consistency 
requirement 7 -C M 2 is fulfilled. 



3. One can arrive at the equations (|3]) without going into the complex p 2 plane. Consider 
only real values of momenta. We look for a point q 2 on the real p 2 -axis such that in the 
neighbourhood of it the function f(p 2 ) can be written in the form f(p 2 ) ~ Zo(p 2 — q 2 + 
ij) + 0(p 2 — q 2 ). The function f(p 2 ) can be expanded around some close-by but arbitrary 
point q 2 , 

f(p 2 ) = Ref(q 2 ) + ilmf(q 2 ) + (p 2 - q 2 ) 



dRef(p 2 ) dlmf(p 2 ) 
+ 1- 



dp 2 



dp 2 



+ 



(5) 
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The ellipses indicate terms of higher order in (p 2 — q 2 ). It follows that 



f(p 2 ) = Z [p 2 -q 2 (l-\) +ll } + ... (6) 
where 7, A and Z are given by 

l(q 2 ) = ^ ~ RfI ' f (7) 
2 RfR'j + If If 

- - m„ 

Thus, if A(g 2 ) is equal to zero one can identify the corresponding value of q 2 with the 
mass of the resonance state and obtain the width of this state from Eq.(^). Indeed, if 
X(q 2 ) = the equality RfR'f + Iflf = is fulfilled. Inserting it in the expression for 7 and 
using that now If R'f — RfI'f = [{R'f) 2 + (If) 2 )If/R'f-> we obtain a result fully coinciding 
with Eqs.(|3]). This shows that when the points q 2 and q 2 — i'j are close to each other an 
approach with the self-consistency condition in the form 

\(q 2 ) = (10) 

is equivalent to the method of the complex pole in the small width approximation. Actu- 
ally, it is a simple consequence of the fact that the derivative of an analytic function f(z) 
at the point z does not depend on the way which the point z has been approached. In 
particular, it does not matter from which direction we calculate the value f'{p 2 ) to obtain 
the result Eqs.@, whether along the real axis or along the imaginary one. The manner in 
which we have arrived at the formula Eq.(jT0|) is not restricted by the requirement 7 -C q 2 . 
It is clear, however, that at large values of 7, (7 ~ q 2 ), the relation of this solution with 
a complex pole location is not so direct. We conclude that in such type of approach, the 
function A plays the crucial role to find a resonance state or to exclude it. 

In the case in which the function f(p 2 ) is real, it is simple to see that A = Rf/{q 2 R'j) 
and the self-consistency requirement leads to Rf — 0. This is a well known fact that the 
mass of physical particles in quantum field theory is associated with the pole of the full 
propagator. 



4. Let us consider now the spectral-function method 0-0. The spectral function is the 
imaginary part of the propagator A(p 2 ) 



£(p) = ImA(p 2 ) = . (11) 

K f + I f 



In this approach it is assumed that the parameters of the resonance can be determined 
from the extremum of the function £(p). The extremum condition is given by 

{R) - I))I) - IRflfR) = . (12) 
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By use of Eqs.(^) and (8) this condition can be written as 




A(? 2 ) + ^¥ =0 • (13) 



It is obvious that this criterium is different from the requirement Eg. (|TTj|) . Therefore, 
one must be fully aware that in the SF approach there is no direct relation with the 
complex pole method and from this point of view can lead to incorrect results. This is in 
fact, very easily understandable. The difference between approaches 2 and 3 and the SF 
approach comes from the fact that in the SF approach one assumes that Z is real. Only 
in that case is the Im A, Eq.(0), maximum at p 2 = M 2 . Indeed if I* — then Eqs. (|T3"D 
and (10) are identical. The origin of a complex Z factor comes from the appearance of a 
background term in the T matrix so that a difference between different approaches signals 
the presence of a non-negligible background. Whenever this happens it is useful to check 
the determination of the resonance parameters via the speed plot technique |J . 



5. Consider now an extended NJL model with spin-1 mesons [II], For the moment, 

let us postpone the rigorous solution of the mass equations and content ourselves to discuss 
the simplest approach to the problem as outlined in section 3. We look for the resonance 
solutions for vector p, axial vector a\ and scalar a mesons, on the basis of Eqs.(|3|). The 
corresponding self-energy functions can be written 

(14) 
(15) 
(16) 





3 

8Gy 


■P 2 J2(P 2 ), 




f al (p 2 ) 


3 

8Gy 


- 6m 2 / 2 (p 2 ) — 


p 2 Mp 2 ), 


Up 2 ) 


m 
AmGs 


+ (4m 2 - p 2 


)Up 2 ) ■ 



Here, Gy and Gs are the constants of the four-quark interaction of the vector, V, and 
the scalar, S, sectors. We also use the symbols m and m for the current and constituent 
quark masses, respectively. The quark-loop integrals I2 and J2 are defined as 



l2 ^ = TT* / ^ / -R{z)e^^\ (17) 

l07T z JO JO Z 

Up 2 ) = ^Y dy(l - y 2 ) [" ^R(z)e^-y 2 \ (18) 
6Lix A Jo Jo z 

The Pauli-Villars regulator R(z) fl2| is chosen in the form 

R(z) = [l - (1 + tzA 2 )e- izA2 ] e~ im2z . (19) 

It corresponds to two covariant subtractions at large momenta p 2 = A 2 ^> m 2 leading to 
a finite result for the p, a\ and a self-energy diagrams. The first two of them contain the 

# 2 Other references can be found in pij. Wc shall follow the notations of that paper here. 
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quadratically divergent integrals. Using the Pauli-Villars method one can avoid these A 2 
terms from the self-energy diagrams. The following properties of the function R(z) are 
important for that: R(0) = R'(0) = 0. 

Evaluating the J 2 and J 2 integrals we get 



16tt 2 



In 



'A 2 + m 2 ~ 



+ 2 



m - m 




(20) 



16vr 2 



1 



+ 



A 2 



1 + 3^ 

2^ VA 2 + m 2 ; v i-l 



Here 



v 



Am 2 ' * 4(A 2 + m 2 ) ' 
and the one-loop function J(x) of the real variable x is equal to 



(21) 



(22) 



J{x) 



In 



1 arctan ■ 



if x < , 
if < x < 1 



(23) 



We remark that our function J 2 (p 2 ) coincides with the corresponding expression from ||, 
precisely we have 167r 2 p 2 J 2 (p 2 ) = 9A 2 Ii(z), where I\{z) is given by Eq.(15) in ||. 

Let us fix the parameters of the model. There are four of them in all: Gs,GV,m 
and A. The constituent quark mass, m, is related to the current quark mass, m, by 
the gap equation of the model. We fix parameters of the model in order to obtain the 
pion decay constant, /„., and the pion mass, m n , close to their physical values, f n ~ 93 
MeV and m v ~ 139 MeV, respectively. The other two input values are determined from 
the vector mesons sector such that one obtains values close to the empirical masses of 
p, m p ~ 770 MeV and a\ meson, m ai ~ 1100 MeV. We investigate two cases, one with 
a large constituent quark mass (m = 390 MeV) leading to a bound state solution for 
the p meson, and one with a small quark mass (m = 180 MeV), where the physical p 
meson appears as a resonant state in the gg-continuum. We obtain for the large quark 
mass case: G<? = 9.41 GeV -2 , G v = 10.27 GeV" 2 , m = 3.9 MeV and A = 1 GeV. With 



these parameters, we find f w = 93.3 MeV, 



139 MeV, 



m r 



776 MeV, m ai 



1042 - z358 MeV, 
0.86 GeV" 2 , G v = 



439 — i0.28 MeV. For the small quark mass case we get: G< 



6.65 GeV 



-2 



m 



lMeV and A = 2.8 GeV. With these parameters, we 



have f n = 93.3MeV, m n = 139 MeV, m p = 674-i300MeV, m ai = 802-i289MeV, m a = 
378 — i2.6 MeV. Note that to determine these resonance properties, we have made use of 
Eqs.(|). 



6. In order to search for resonances which lie in the complex plane on the second Riemann 
sheet, we introduce the following analytic continuation of the functions J(£) and </(£). 
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The function J(£) has a branch point Pq at z = 1 and a cut along the real axis from P 
to +00. Analytically continuing this function on the second Riemann sheet we obtain the 
following analytic function J-'oi 2 )^ where z G C, 



Mz) 




-m + log 



1 + t/l 



1 - a/1 



(24) 



We are looking for resonances whose masses are smaller than the cut-off A so that one 
does not cross the cut of the function J(£) which lies from Pi to +00 where Pi is at z = 1, 
that is at z = 1 + (A 2 /m 2 ). This function is thus given in the complex plane C by the 
analytic function T\(z) 



where z = zm 2 / (A 2 + m 2 ). Let us consider now the mass equations 

f a (z) = 0, a = <7,ai,p . (26) 
It is knownp 3 ! that in the case of cr or ai mesons there exists only the resonance solution 



of these equations. To show them we shall plot functions u a (Gy) and in some cases also 
r) a (Gy), where (uj a — ir] a ) 2 = 4m 2 z a and z a is a solution of the corresponding equation in 
Eq.(^) at some value Gy. We shall call these functions a resonance branch of Eq. (|26|) if 
they are a mapping of resonance solutions of this equation at different values of Gy. The 
corresponding formulae are 



mK aj rj a = -2mK a l lmz a , K a = jRez a + \z a \. (27) 



By our definition z is proportional to the complex p 2 . In turn complex p 2 reads p 2 = 
(po — ii]'/2) 2 — p 2 . In the rest frame p = and then p 2 = (M — irj/2) 2 . Thus u has a 
meaning of mass and 77 is related to the width of the resonance state. If 77 is small one 
has p 2 = M 2 — ij, where 7 = Mr/. We used this fact in paragraph 2. 

We plot resonance branches for p and a\ mesons in Figs. 1-2 for two different values 
of quark masses: m = 390 MeV and m = 180 MeV, using the same set of parameters 
(m,m n ,A) as in section 5 and varying Gy. Notice that in the case of the p meson only 
the bound state branch has been known until now. At p 2 > Am 2 this curve smoothly 
transforms into a "virtual state" section. This has been first observed in ||. In nonrel- 
ativistic scattering theory, the appearance of a "virtual state" is directly related to the 
absence of a near-threshold resonance in the s-wave. In the relativistic theory, as we have 
found here by studying the NJL model, the resonance branch appears together with a 
bound state- "virtual state" branch as a separate curve. It is never smoothly connected 
with the bound state branch and it has rather large values of the imaginary part r\ p . For 
large values of quark mass m ~ 400 MeV (see Fig.l) the p meson solution belongs to the 

# 3 One can also prove this fact analytically. 
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bound state branch. However, at relatively small values, m ~ 200 MeV, as one can see 
from Fig. 2, the p meson solution is a resonance state in the quark-antiquark continuum. 

We plot in Figs. 1-2 only the real parts, u a , of resonance solutions m a = u a — ir) a . 
Therefore the intersection of the resonance branch with the "virtual state" in Figs. 1-2 
is only a result of this projection on the real m a axis. Really the crossing curves have 
different imaginary parts r] a in these points. We draw them in Figs. 3-4. For m = 390 
MeV we obtain for the empirical value of f n = 93.3 MeV the p-meson bound state mass 
m p = 776 MeV and the complex mass for a\, m ai = 1000 — i263 MeV. For m = 180 
MeV, the complex masses of resonance states which correspond to the empirical value of 
U = 93.3 MeV are m p = 674 - i292 MeV, and m ai = 802 - i289 MeV. These values are 
a bit smaller than the experimental ones. To clarify the dependence of the result we 
indicate in Figs. 1-2 boxes inside which the f n value changes from f n = 100 MeV on the 
left side to = 92 MeV on their right side. The physical solutions lie inside this box. In 
the other figures we directly plot f v as a function of Gy. As discussed in section 4, we 
have also checked our results via the speed plot technique. Very good agreement between 
the two methods was obtained. 

It is also interesting to compare the results of the exact solutions of the complex BS 
equation for the p and the a\ with the approximate estimations of section 5 (see Figs. 3-4). 
We start the comparison for the case of a large quark mass m = 390 MeV in Fig.3. As 
one can see the expansion approach works quite well for the a\ meson. The deviation in 
the results of the two approaches is smaller at higher values of Gy where rj ai is smaller 
as it should be expected. We also plot here the r\ p values. It follows from section 3 that 
at p 2 < 4m 2 , Eqs.(^) describe a bound state branch. At large values of the quark mass, 
we have Re(z p ) <C 1, far away from the region of applicability of Eqs.(|3]) to resonance 
solutions. The resonance branch as an approximate solution cannot be found in this case. 
The physical branch for p is a bound state solution, for which the exact solution fully 
coincides with the solution of Eqs. ([$]). In Figs. 4a and 4b we compare the full results 
with the expansion approach at small quark mass values m = 180 MeV in the case of the 
ai and p mesons respectively. In the case of the p-meson the larger deviations observed 
between approximate and exact solutions at Gy > 12 GeV -2 values stem from the fact 
that Re(z p ) < 1. In this region, 12 < Gy < 20 GeV -2 , the expansion curves interpolate 
between the resonance solutions and the bound state branch (if it exists). In Fig.4b we do 
not plot the bound state branch which starts at Gy = 20 GeV -2 . As mentioned above, 
Eqs.(|3|) are not appropriate to describe resonance solutions in this region. The exact 
scalar resonances coincide with the ones obtained with the approximate ones (section 5) 
and are practically constant as functions of Gy. 

7. As in ref. M we have studied here the resonance solutions in the extended NJL model 
in the quark-antiquark continuum. There is only one difference in our Lagrangian in 
comparison with that paper. We consider the model without the 't Hooft determinantel 
interaction. This interaction is important for solving the U(l) anomaly problem and, 
therefore, essential for the calculation of r] — rj' meson masses |R|. Since we restrict our 
consideration to the case of p, a\ and a mesons one can neglect these six quark vertices. 
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In any case, the results found in our study are independent of the precise form of the 
underlying NJL-type model. In ref.||, the authors could not find the resonance solution 
in the case of the vector p meson. Instead, the "virtual state" has been found. In this 
letter we have studied this question more thoroughly. We have shown that the relativistic 
BS equation for the p meson state possesses also a resonance solution. This solution 
has a rather large imaginary part and corresponds to the physical p at small values of 
constituent quarks masses m ~ 200 MeV. The imaginary part is a consequence of the 
deconfinement property of the NJL model and should be considered an artefact of the 
model. Nevertheless, the existence of resonance solutions is a very important feature of the 
model. It means in particular that the p meson state in the model exists as an asymptotic 
field and contributes at large time even at small values of m, although the probability 
for the detection of the particle decreases exponentially with a characteristic decay time 
r(p) = uj(p)/'~f. The well-known Lagrangian (bosonization) approaches to the NJL model 
|| use this fact without proof. An asymptotic p meson field is usually introduced in the 
bosonized NJL Lagrangian by the replacements of variables in the generating functional 
and finally, after the fixing of parameters, it happens often that the quark mass value is 
of the order of 200 -j- 300 MeV. This procedure has a justification only if this state really 
exists in the model as a resonance in a quark-antiquark continuum. Our work proves for 
the first time this important fact. 
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Figure captions 



Fig.l Dependence of the real part u a , (a = p, and ax) of the exact vector, p, and axial- 
vector, ai, resonance solutions on Gy at large values of constituent quark masses 
m = 390 MeV. The solid line corresponds to the bound state branch. It smoothly 
touches the qq threshold line, 2m. The dashed lines indicate resonance branches of 
p and di mesons. The dotted line is a "virtual-state" solution. The physical p is a 
bound state solution. The solid box indicates the region bounded by the values of a 
pion decay constant between /„. = 100 MeV (left side of the box) and /„- = 92 MeV 
(right side) where the physical p meson lies. 



Fig. 2 Same as Fig.l for the small constituent quark mass m = 180 MeV. 



Fig. 3 Comparison of the exact solutions (dashed lines) of the a>\ meson BS equation for 
the quark mass m = 390 MeV at different values of Gy with the corresponding 
solutions of the Eqs.(|) (dotted lines). Indicated is also the exact imaginary part r\ p 
and /„■. For more details see text. 



Fig. 4 (a) Comparison of the exact solutions of the a\ meson BS equation for the quark 
mass m = 180 MeV as functions of of Gy with the corresponding solutions of the 
Eqs.(|3|), same notation as in Fig. 3. (b) Same as (a) for the p meson. Indicated is 
also f n . 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4a 
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Figure 4b 
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